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We investigate, within the fluctuation-exchange approximation, a correlated-electron model for 
Y2Ba4Cu70i5 represented by two inequivalent Hubbard layers coupled by an interlayer hopping 
t±. An energy offset 5 is introduced in order to produce a different charge carrier concentration 
in the two layers. We compare several single-particle and magnetic excitations, namely, the single 
particle scattering rate, the spectral function and the spin lattice as well as spin-spin relaxation 
times in the two layers as a function of 5. We show that the induced interlayer magnetic coupling 
produces a tendency to "equalization" of the magnetic properties in the two layers whereby antifer- 
romagnetic fluctuations are suppressed in the less doped layer and enhanced in the heavily doped 
one. The strong antiferromagnetic bilayer coupling causes the charge carriers in the plane with 
larger doping concentration to behave similar to those of the underdoped layer, they are coupled to. 
This effect grows for decreasing temperature. For high temperatures or if both layers are optimally 
or overdoped, i.e. when the antiferromagnetic correlation length becomes of the order or smaller 
than one lattice site the charge carrier and magnetic dynamics of the two layers is disconnected 
and the equalization effect disappears. These results are in good agreement with NMR experiments 
on Y2Ba4Cu70i5 by Stern et al. Phys. Rev B 51, 15478 (1995). We also compare the results with 
calculations on bilayer systems with equivalent layers as models for the constituent compounds 
YBaaCugOr and YBa2Cu408. 



PACS numbers; 74.72.-h,71.27.-|-a,76.60.-k 



I. INTRODUCTION 

There is a large amount of consensus that the anoma- 
lous properties of cuprate superconductors are caused 
by the strongly interacting electrons within the Cu02 
planes. However, in particular the observation that the 
highest superconducting transition temperatures belong 
to compounds with more than one layer per unit cell ini- 
tiated various investigations of out-of-plane properties. 
The observation of a rather strong coupling between ad- 
jaceut layers has been made by inelastic jaeutron scatter- 
ing^l (INS), nuclear magnetic resonanceoH (NlVffi) and 
indirectly also in Raman scattering experimentsO. Fur- 
thermore, the observation of a qualitatively different be- 
havior of the odd and even channel in INS including a. 
sharp resonance feature, found solely for odd excitationsLJ 
and of a bilayer splitting of the Fermi surface foundJ: 
angular resolved photoemission experiments (ARPES)l3' 
demonstrate that low energy excitations of cuprates are 
affected by the presence of more than one layer per unit 
cell. Related to these issues is the interesting question 
of the c-axis transport and, the, occurrence of a c-axis 
Josephson plasma excitationld'tj, which may turn out to 
be a new probe of the vortex statics and dynamics of the 
superconducting state. 

A very interesting perspective on the nature of the cou- 
pling between Cu02-layers was offered by NMR experi- 
ments by Stern et al. on Y2Ba4Cu70i5 (247). This ma- 
terial has a variety of structural similarities to the exten- 
sively studied YBa2Cu307 (123) and YBa2Cu408 (124) 



systems. The main difference in the crystallographic 
structure of 123 and 124 is the double CuO chain in the 
latter. The compound 247 can be considered as a nat- 
ural multilattice, consisting of alternating 124 and 123 
blocks. The bilayers in 247 are correspondingly build up 
of one Cu02 layer which belongs to the 123 block and 
one layer of the 124 block. Based on the analysis of the 
NQR spectra it turned out that the charge carrier con- 
tent in these nonequivalent adjacent layers is very close to 
that of the related parent compounds of the two blocks, 
i.e. one plane has a similar charge carrier concentration 
to the slightly overdoped 123 system whereas the other 
layer corresponds to the underdoped 124 system. Inter- 
estingly, the highest transition temperature (Tc = 95 K) 
is for 247, which has to be compared with the 92 K for 
123 and 82 K of the 124 system. The main experimental 
observations of Ref. ||,^ are the following: (i) the low- 
temperature Knight-shift suppression for both planes is, 
despite their different charge carrier concentration, simi- 
lar and behaves like in the underdoped 124 system, even 
though the high temperature values of the Knight shift 
are the same as in the corresponding 123 and 124 com- 
pounds, (ii) both ^^Cu spin lattice relaxation times of the 
two different planar Cu-sites show a spin pseudogap in 
l/^^TiT, even though it is barely present in the 123 par- 
ent compound itself, (iii) the interplane transverse relax- 
ation rate, as measured in a spin-echo double resonance 
experiment (which characterizes the interplane magnetic 
susceptibility) increases for decreasing temperature faster 
than the intraplane relaxation rate. Thus, the main con- 
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elusions from these observations are that, for high tem- 
peratures, the two planes are rather disconnected and 
behave similarly to their parent compounds, whereas for 
decreasing temperatures, the increasing interlayer mag- 
netic coupling enforces even the slightly overdoped plane 
to behave like an underdoped system. 

For a proper interpretation of these interesting exper- 
imental data and, in a more broader context, for a bet- 
ter understanding of the bilayer coupling in cuprate su- 
perconductors in general, it is essential to investigate to 
what extent one can describe the main trends of these 
data within a model of coupled layers, only different by 
their charge carrier concentration or whether one needs 
to make qualitatively new assumptions about the nature 
of the bilayer coupling. 

One promising approach for the description of bilayer 
phenomena is based on a Hubbard Hamiltonian with lo- 
cal repulsive Coulomb interaction, where the interplanar 
coupling is caused solely by an interplane hopping ele- 
ment t±. This model, restricted to the case of equiva- 
lent layerS|-(^f^ 0), has been investigated within various 
techniqueal3iiJ. Additional insight can be gained using a 
self-consistent summation of bubble and ladder diagrams 
(fluctuation exchange aBjpoximation). The main results 
of these investigational3il3 are enhanced antiferromag- 
netic spin fluctuations due to layer coupling causing, in 
bilayer systems as well, a (1^2 _y2 symmetry of the super- 
conducting order parameter, a predominantly incoherent 
low energy c-axis charge transport even though the bi- 
layer splitting stays intact, and an enhancement of the 
relative strength of interlayer vs. intralayer coupling for 
decreasing doping. 

In this paper, we additionally consider the effect of an 
energy offset 5, which produces a different charge carrier 
concentrE^tipn in the two Hubbard layers (cf. also Scalet- 
tar et al.t3). This is a suitable model to describe the 
peculiarity of the 247 compound Y2Ba4Cu70i5, whose 
bilayers are build up of one Cu02 layer belonging to 
an YBa2Cu307 block and one to a YBa2Cu408 block. 
We evaluate several single-particle and magnetic exci- 
tations, namely, the single-particle scattering rate, the 
spectral function and the spin lattice as well as spin-spin 
relaxation times in the two inequivalent layers as a func- 
tion of 6 within the fluctuation-exchange approximation. 
We show that the interlayer coupling produces a ten- 
dency to equalization of the antiferromagnetic properties 
in the two layers whereby antiferromagnetic fluctuations 
are suppressed in the less doped layer and enhanced in 
the heavily doped one. This equalization effect turns out 
to be enhanced in the presence of antiferromagnetic fluc- 
tuations in the system and to be almost absent when 
the antiferromagnetic correlation length becomes of the 
order or smaller than one lattice site and to ultimately 
decrease for increasing temperature. These results are in 
good qualitative agreement with NMR experiments on 
Y2Ba4Cu70i5 by Stern et aluu. We also compare the 
results with calculations on bilayer systems with equiv- 
alent layers as models for the constituent compounds 



YBa2Cu307 and YBa2Cu408. 

A first theoretical investigation of the experimental 
findings of Ref.B has been given by Millis and MonienE3, 
who could determine the size of the interlayer exchange 
coupling from an analysis of the interlayer cross relax- 
ation time. These authors also discuss that the 41 meV 
excitation observed in superconducting YBa2Cu307 is a 
collective mode pulled down below the superconducting 
gap by interactions, and that the observed antisymme- 
try under interchange of planes follows from the non- 
negligible value of J±. An analysis of the coupling be- 
tween an undoped layer and an underdoped one, simi- 
lar inpSpirit to ours, has been carried out by Scalettar 
et al.ta. These authors study the pairing mechanism, 
which arises from the coupling of holes in doped layers 
to spin fluctuations in the undoped layers in analogy with 
the Ginzburg-type scenario for the coupling of electrons 
through excitons in a doped semiconductor. However, it 
turns out that magnetic fluctuations in the undoped layer 
are strongly suppressed by the coupling with the doped 
layer and superconducting correlations are reduced by 
the interplane coupling at least at the temperatures ac- 
cessible to the simulations. The study of the coupling 
between a strongly antiferromagnetic and a doped sub- 
system has some similarities with the "stripe scenario" 
where hole-poor antiferromagnetic regions are considered 
to be in contact with hole-rich superconducting regions. 

The paper is organized as follows: In Sec. || we present 
our model for coupled layers with different charge car- 
rier concentration and summarize the main concept of 
the fluctuation exchange approximation, used for the ap- 
proximate investigation of the model. In Sec. HI we 



present our numerical results with particular emphasis to 
the single-particle and magnetic fluctuations in the two 
layers and focus on the anisotropy and on the tendency 
of equalization of this effects. In order to make contact 
with the experimental investigations on the 247 system, 
we discuss at length the temperature dependence of var- 
ious NMR quantities in Sec. IV. Finally our results are 
summarized in Sec. 



II. MODEL AND TECHNIQUE 

In order to describe the strong electronic correlations 
in the high-Tc superconductors and the particularities of 
the system Y2Ba4Cu70i5, consisting of two layers with 
different charge carrier concentration, we use a system 
of two two-dimensional Hubbard layers coupled by an 
hopping element t± . After Fourier transformation of the 
intraplane sites into momentum space with in-planc mo- 
mentum k, the Hamiltonian reads: 



k,o- 



where cj^ ^ creates a particle with spin a at momentum 
k in layer li. Furthermore, n^^i is the density operator 
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at lattice site and spin t, and -ffo(k) the Hamilton 
matrix for the noninteracting system 

In order to describe theoretically a different charge carrier 
concentration in the two layers, we additionaly introduce 
an on-site energy 6 in the second layer, effectively modify- 
ing its chemical potential. The planes are coupled solely 
through a bare interplane hopping t±. Furthermore, the 
bare energy dispersion in each plane is 

eic = — 2i (cos kx + cos ky) — At' cos kx cos ky 

~2t" [cos(2A:^) + cos(2fcj^)] . (3) 

thus including second and third-neighbor hopping pro- 
cesses {t',t") to better model the Fermi surface for the 
system under consideration. In theJaUpwing calculation, 
we always set tj_/t — 0.4, U — 4mM3 and measure the 
energies in units of the next-nearest neighbor hopping t. 

A diagonalization of ffo(k) leads to the bonding and 
anti-bonding bands of the noninteracting system 

e±^.k + ^±y^+ii. (4) 

The single-particle excitations and the thermodynamic 
properties of the interacting system are deduced from 
the Green's function G(k, iajm) obtained through Dyson's 
equation which for a two-layer system generalizes to a 
(2 X 2) matrix equation 

G"^(k,ia;„i) = (icjm + - i?o(k) -^(k.iuj^). (5) 

Approximations are introduced by the explicit choice of 
the self-energy S(k, zcj,n). Here, we use the expression for 
the self-energy given by the FLEX approximationEll with- 
out particle-particle vertex contributions. Within the 
FLEX, the irreducible particle-particle vertex is solely 
the repulsive Coulomb interaction U and consequently 
irrelevant. Interference effects between the particle- 
particle and particle-hole channel, which may be of rel- 
evance for an understanding of the pseudogap state of 
underdoped cuprates at low temperatures, are beyond 
the scope of this paper. 

Introducing the shorthand notation k = (k, zwni) and 
q = (q, zj/n) for convenience, the matrix for the self- 
energy E/;/ (fc) reads 

^wik) = ^Y.^iAk-k')Gu,ik'), (6) 

where (3 — -j^^ is the inverse temperature, and the ef- 
fective interaction Vw results from an infinite series over 
spin- and charge-fluctuations and is given by 

V{q) = —{l-Ux{q))-\{q) 

+ ^il + Uxiq))-'xiq) (7) 
-U\{q) . 



Note, that V and x are (2 x 2) matrices, i.e. matrix in- 
version and multiplications have to be used. The bare 
particle-hole bubble xii'{q) consists of dressed Green's 
functions 

Xw (9) = - ^ E G"' + 9)^''' (^^) ■ (8) 

In Eq. |[ the Green's functions are determined self- 
consistently by solving the set of coupled equations Eqs. 
(^,||j7|j|). During the self-consistency cycle, we fix the on- 
site energy 5 and the particle number rii — \ — xi oi the 
first layer, while the chemical potential /i and the particle 
number of the second plane n2 — \ — X2 are determined 
at each step. It turns out that the total particle number 
n = ni + does not essentially depend on temperature, 
which makes the physical interpretation of our numerical 
results more straightforward. To avoid the uncertainties 
related to a numerical analytical continuation of corre- 
lation functions from the imaginary Matsubara to real 
frequencies, we use the recently proposed real- frequency 
approach to the FLEX approximation.E3 

In the following, we shall mainly focus our attention on 
two different parameter sets in order to mimic a situation 
with strong and weak antiferromagnetic fluctuations, re- 
spectively. Specifically, we use a parameter set, for sim- 
plicity labeled by "A" , with t' = -0.38i, t" = -0.06i and 
a second one, labeled by "B" , with t' = -0.20i, t' = O.lSt. 
As will be shown below, the parameter set A corresponds 
to a system with pronounced antiferromagnetic fluctua- 
tions at low temperature, while B has much weaker ones. 



III. RESULTS 



In order to investigate the effects of different charge 
carrier concentrations in coupled bilayer systems, we 
start the discussion with the 5 dependence of the doping 
of the second plane, X2 {5) . Here, the doping of the first 
plane xi and the on-site energy 5 are independent vari- 
ables for us, while X2 comes out from the self-consistent 
calculation. In Fig. ^we present X2{5) for both parameter 
sets A and B for different doping levels of the first plane 
as a function of 5. From Fig. ^ it can be seen that the 
hole doping of the second plane is directly proportional 
to its on-site energy 5. Note that the constant of propor- 
tionality, i.e. the slope dx2/d5^ is almost independent 
on t' ,t" and xi, because the particle number difference 
5n — ni — n2 between both planes is mainly governed 
by the energy dependence of the effective "chemical po- 
tentials" ^ and fj, — S. On the other hand, correlation 
effects do play a role here, since the slope depends on 
U/t. Note that this energy difference is determined self- 
consistently in our theory and therefore does not only 
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depend on 5, but varies, due to self-energy renormaliza- 
tions, also with the strength of the Coulomb interaction 
U. The reason why we follow the strategy of keeping the 
doping of the first plane fixed and only change the dop- 
ing of the second plane is that we want to investigate a 
possible induced coupling effect between the planes. The 
question is whether single-particle and magnetic proper- 
ties of the first plane are influenced by the doping of the 
second plane. 
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FIG. 1. Hole density of the second layer as a function of 
its on-site energy S for parameter sets A and B (described in 
the text) and different dopings of the first layer [T = 0.02t). 
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FIG. 2. Real part of the diagonal elements 

of the self-energy, ReE;i(k,tj) for k = (7r,0) and 
xi = 0.08, T — 0.02i. (a) and (b) correspond to the parameter 
set A and (c) and (d) to the set B, respectively. 

In Fig. we show the real part of the diagonal el- 
ements of the self-energy Y,ij(k,uj) with k = (tt, 0) for 
both layers and both parameter sets A and B. The dif- 
ferent curves in each panel are for various on-site energies 
S. Fig. ^(a) displays ReS for the first (less doped) layer 
with parameter set A, i.e. for t' = —0.38t, t" = 0.06t. 
This figure shows that ReSn indeed depends on S al- 
though the hole concentration is not changed by 6 in the 
first layer. This demonstrates that feedback effects due to 
the interlayer coupling modify also the properties of the 



plane where the charge carrier concentration is kept con- 
stant. The dip- like structure for small 5 is a precursor of 
new quasiparticle states on the shadow of the Fermi sur- 
face due to strong antiferromagnetic fluctuations. This 
effect, which is strongest for the case of equivalent planes 
(S — 0), has been discussed in Refs. [l7|Jlq . An increase 
of (5 and hence of the total hole doping of the system 
leads to a decrease of the interplanar antiferromagnetic 
coupling and thus of the dip structure. Panel (b) in 
this figure shows ReS22((7r, 0), for the second (heav- 
ily doped) layer. Since it is this plane which is primary 
altered by 8 the influence of changing the on-site energy 
(5 is, as expected, considerably more pronounced. These 
results show that the two layers are strongly connected 
and a change in carrier concentration of the second layer 
strongly influences the single-particle properties of the 
first layer as well, although the doping is unchanged here. 
However, whether the two planes are connected or not de- 
pends on the values of the parameters of the model. For 
example, a completely different situation is found for the 
parameter set labeled by B. For this choice, a variation 
of 8 influences the second layer [panel (d) of Fig. |^ , but 
has no effects on the self-energy of the quasi-particles in 
the flrst layer [panel (c)], indicating independent planes. 
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FIG. 3. Imaginary part of the diagonal elements 
of the self-energy, ImEu(k,aj) for k = (7r,0) and 
xi = 0.08, T = 0.02t 



The connection between the two planes is also visible 
in the scattering rates which are related to the imaginary 
part of the self-energy and presented for both layers and 
both parameter sets in Fig. ^. For the imaginary part of 
the self-energy, which is more sensitive to low energy ex- 
citations, we observe an even closer connection between 
the two planes for parameter set A than for the real parts. 
In addition, these figures demonstrate that precursors of 
a spin density wave state around (tt, 0), for low (5, are 
rather incoherent due to the strong scattering rates at 
these energies. Note also, that even the quasiparticles at 
the chemical potential (cj = 0) suffer strong scattering, as 
indicated by the rather large values ImI]H((7r, 0), w = 0). 
Similarly to the real part, the changes in the scatter- 
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ing rates caused by 5 are rather moderate for parameter 
set B, even in the second layer which is directly altered 
by 6 through its doping. Thus, both the real and the 
imaginary parts of the self-energy suggest a tendency of 
equalization between the layers for parameter set A, but 
not for B. 

Nevertheless, the two planes turn out to be discon- 
nected for regions far from the Fermi surface, e.g. close 
to k = (0, 0), even for data set A. 
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FIG. 4. Real and imaginary part of the two diagonal ele- 
ments of the self-energy for parameter set A and for k — (0, 0). 

The real and imaginary parts of the self-energy for pa- 
rameter set A and k = (0, 0) are presented in Fig. |^. As 
can be seen in part (a) and (c) of this figure, changes 
in the on-site energy S of the second plane have almost 
no effect on the first plane. The two planes are thus con- 
nected only for momenta close to the Fermi surface which 
are strongly affected by antiferromagnetic fluctuations at 
low energy and in particular for hot quasiparticles states 
around (tt, 0). 
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FIG. 5. Scattering rate in 

the first plane — ImEii(k, c<j = 0) along the standard path 
in the Brillouin zone for parameter set A with xi = 0.08, 
T = 0.02t, S/t = 0.0, 0.2, 0.4. The arrows indicate crossing of 
the anti-bonding (kp") and bonding (kj? ) Fermi surfaces (see 
also Fig. 0). 

To further elucidate the momentum-resolved equaliza- 



tion effects for parameter set A, we show in Fig. ^ the 
scattering rates at the Fermi energy for the first plane, 
namely — Iml]ii(k, w = 0) along the standard path in the 
Brillouin zone, for different values of the on-site energy 
S of the second plane. This figure demonstrates that 
the scattering rates are strongly modified at the Fermi 
surface with strong effects in the regions close to (7r,0). 
On these regions, the large number of states associated 
with the flat bands produces strong scattering processes, 
whenever the interaction connects the van Hove regions 
of the bonding and antibondig band. We stress again 
that the variation of the first plane is solely caused by 
its correlation with the second plane since the doping of 
the first plane is kept constant. The corresponding real 
part of the self-energy (not shown) also reveals coupling 
effects for all k close to ki;-. 
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FIG. 6. Effective interplane hopping t']_ (k, uj) as a func- 
tion of the onsite energy 5 for t± = OAt, xi = 0.08, T = 0.02t 
for parameter set A. The inset shows the results for parameter 
set B on the same scale. 
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FIG. 7. Effective interplane hopping t'f*(k, w) as a func- 
tion of temperature for parameter set A with xi = 0.08, 
S — 0.0. Note the stronger dependence on T than on 5 when 
compared to Fig. ^ 



We now address the question of how much the hop- 
ping tj_ between the planes is affected by the interaction. 
Dyson's equation (Eq. ^ suggests a momentum- and 
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energy-dependent effective interlayer hopping according 
to 



tf{k,uj) =ti +ReEi2(k,w) . 



(9) 



Like the real and imaginary parts of the diagonal ele- 
ments of I](k, w), this quantity is strongly affected for 
k « (tTjO), therefore, we restrict ourselves to this mo- 
mentum. 

In Fig. ^ we show t°^{'k,uj) for parameter set A at a 
fixed temperature T — 0.02i for a series of on-site ener- 
gies 6, while the inset presents the same quantities for 
parameter set B. Looking at Fig. ^, we observe that t± 
is strongly renormalized close to the chemical potential 
(w = 0), and even more important, that the antiferro- 
magnetic fluctuations again enhance the interlayer con- 
nection, in this case represented by the effective hopping 
. On the contrary, t± is essentially unrenormalized, 
independently of the value of 6, in the presence of the 
weak antiferromagnetic fluctuations in parameter set B 
(see inset of the figure) . This result seems to stress once 
more the fact that the planes are strongly connected for 
the "antiferromagnetic" parameter set A, whereas they 
are essentially independent for the parameter set B. We 
thus focus our attention on parameter set A where the 
eff'ects are stronger. For (5 = the effective hopping 
i'^((7r, 0), 0) « 0.8t is roughly twice as large as the bare 
t± — OAt. Thus, the hopping between the planes is am- 
plified rather than blocked by electronic correlations. 

Fig. 1^ also shows that the renormalization of tj_ de- 
creases with increasing on-site energy d, which causes the 
second plane to be less magnetic. However, an increase 
of 6 has a surprisingly weak effect on t]^((7r, 0), ui) when 
compared with an increase of the temperature T, which 
is shown in Fig. |7[ This strong temperature dependence 
of ^^^((Tr, 0),a;) shown in the figure suggests that the in- 
terlayer coupling is related to a small energy scale ui*. 
For high temperatures, the thermal fluctuations destroy 
the correlations on the small energy scale uj* . 
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FIG. 8. Spectral functions ^(k,tj) of the bonding and 
anti-bonding bands at k = (vr, 0) for parameter set A with 
xi = 0.08,5 = 0,T = 0.02t 



Although the above discussion suggests an enhance- 
ment of t± due to interaction effects, the band split- 
ting between the bonding and the anti-bonding band 
A(k) = uj^ — uj^ goes in the opposite direction and 
is reduced with respect to its bare pJufii Aq = 2ij^, in 
agreement with previous conjecture^JtJ. Fig. ^ shows 
the spectral functions A(k, uj) of the bonding and anti- 
bonding band for k — (tt, 0) and indicates a renormal- 
ized band splitting A(7r, 0) « OAt = 0.5Ao. Thus, while 
the interlayer hopping seems to be enhanced by about a 
factor of 2 at low energies in the presence of strong an- 
tiferromagnetic fluctuations, the band splitting behaves 
in the opposite way and it is reduced by about the same 
factor for this parameter set. This different behavior be- 
tween and A(k) is rather surprising, although it may 
be understood by the following argument. On the one 
hand, the quasiparticle interplane hopping without resid- 
ual interaction, related to the off-diagonal energy term 



, is enhanced due to the fact that quasiparticle of two 



neighboring sites on the two planes are nearly antiferro- 
magnetically ordered and thus have a larger amplitude to 
hop. On the other hand, the whole hopping amplitude, 
related to A(k), is suppressed (by a larger factor than 
the enhancement of f^^), due to the Hubbard repulsion 
U. 

Equalization effects between the planes are also ob- 
served in two-particle quantities like the spin response as 
deduced from the spin-spin correlation function. This is 
given (in the layer representation) by 



X 



= (q,c^) = 2[l-C/x(q,c^)]"'x(q,^) 



(10) 



The static spin-spin correlation function xf;^(q, = 0) 
along the standard path (0, 0) — > (tt, 0) (tt, tt) (0, 0) 
in the Brillouin zone is shown in Fig. |^ for I — 1 and for 
I — 2 for parameter set A and B. 
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FIG. 9. Static spin susceptibility Xiti^l,^ ~ 0) along the 
standard path in the Brillouin zone for the two layers and 
both parameter sets A and B. (T = 0.02t, xi = 0.08) 

In the case of parameter set A, the spin response in 
both planes is strongly peaked at q = (tt, tt) indicating 



6 



considerable antiferromagnetism in the Hubbard planes. 
However, even more important is the strong dependence 
of Xii (q, 0) on S which is solely due to the interplane cou- 
pling since the doping in the first plane is fixed. In agree- 
ment with the effects observable in the single-particle 
spectrum represented by the self-energy, this clearly re- 
veals a strong connection between the planes. Fig. ||(a) 
and (b) also show that the spin response is very sensi- 
tive to a variation of 6 and the antiferromagnetism is 
suppressed if S is increased. This is due to the increasing 
hole concentration in the total system which tears it away 
from half-filling where antiferromagnetism is strongest. 
The static spin response for the data set labeled by B is 
shown in Fig. ^(c) and (d). For this choice of parameters, 
the first plane is again completely disconnected from the 
second one since a change of 6 only influences the spin 
response of the second plane while that of the first plane 
is not affected at all. A clue towards the understand- 
ing of the difference between the parameter sets A and 
B is already found in the behavior of the spin response. 
Comparing Fig. ^ (a)-(b) with Fig. || (c)-(d) reveals a con- 
siderable (one order of magnitude) smaller value for the 
static spin response at (tt, tt) in the second case. Further- 
more, the spin correlation length ^, which is the inverse 
of the half width at half maximum of x^^il^ 0) is much 
smaller (of the order of 1 lattice spacing) for parameter 
set B. This implies, as expected, that the interplane con- 
nection is intimately related to strong antiferromagnetic 
correlations in the planes. 
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FIG. 10. Bonding and anti-bonding Fermi surfaces for 
parameter sets A (left panel) and B (right panel) with 
= 0.08,5 = 0.2,r = 0.02t. 

The strength of the antiferromagnetic fluctuations are 
in turn quite sensitive to the shape of the Fermi surface, 
especially if large regions of the Fermi surface can rbe 
linked by the antiferromagnetic momentum Q = (tt, 7r)E3. 
Keeping this in mind, the differences between parameter 
set A and B are caused by their different Fermi surfaces. 
These are shown in Fig. |l^. While the antibonding and 
bonding Fermi surfaces associated with parameter set 
A are closed around (0,0) and (7r,7r), respectively, both 
Fermi surfaces of B are closed around (tt, tt). Even more 
important is the fact that large regions of both Fermi sur- 
faces in case A may be connected by Q. These regions 
are close to (tt, 0), i.e. close to the van Hove singularities 



in the density of states, thus opening various channels for 
antiferromagnetic scattering processes. 



IV. RELATION TO EXPERIMENTS 

We now turn to the question how our theoretical 
calculations compare with experimental results. We 
thus concentrate our attention to NMR experiments on 
Y2Ba4Cu70i5 performed by Stern et alBB. Here, the 
experimentally relevant quantities are the spin-lattice re- 
laxation time Ti and the Gaussian contribution T2G to 
the nuclear spin-spin relaxation time T2. _ 

As pointed out by ShastryEZi and Mila and RiceES, the 
spin-lattice relaxation time Ti is related to the spin sus- 
ceptibility x^^, via the expression: 



TiT ^^0 2N h ^ hu 
q 



(11) 



where fc(q) is the form factor resulting from the Fourier 
transform of the hyperfine interaction 



Fc(q) = {Aab + 2B[cosq,r. + cos qy]} 



(12) 



Thus, Ti probes the slope of the imaginary part of 
X^^(q, w) for w — » 0. In contrast to this, the Gaussian 
component of the transverse relaxation time T2G depends 
on the static susceptibility and is given by 



rp — 2 

^2G 



0.69 



l^i^2^(q)x^^^(q,0) 



l5]F,ff(q)x^^(q,0) 



(13) 



as pointecLout independently by TakigawaO and Thelen 
and Pineg23. While Eq. |l3| applies for the diagonal ele- 
ments of Xif 7 i-6- for the in-plane relaxation rates, the 
corresponding inter-layer relaxation rate 1 /T2Q needs not 
to be corrected by the_self-interacting hyperfine interac- 
tion and thus reads :B'E2l 



2GJ 



0.69 1 ^ 
q 



^eff(q)[xf2(q,o)]' (14) 



The form factor i^eff(q) in the last two equations is sim- 
ply obtained form Eq. (^ by replacing Aab with Ac- Al- 
though we are studying inequivalent Hubbard planes, we 
assume for simplicity that the hyperfine constants, which 
are usually extracted from Knight shift experiments, are 
identical in both planes. This is supported by the rather 
moderate variations of the hyperfine coupling constants 
for different cuprate superconductors. Quantitative dif- 
ferences may be obtained if one takes into account differ- 
ent hyperfine constant, although the temperature depen- 
dence should not change. To be specific, we adopt here 
the values recently given in an analysis of NMR expqp-, 
iments on YBCO and LSCO by Barzykin and PincsEil 



7 



and set Aab = 0. SAB, Ac = —AB and the energy scale 
B = 3.82 X 10"'^ eV. Note, that since NMR probes the 
local environment of the spins, all momenta q contribute 
to the relaxation times although the main contributions 
come from the regions q ^ (ti", tt). Hence, the behavior 
of the NMR relaxation times are strongly influenced by 
the antifcrromagnetic response of the system. 
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FIG. 11. Temperature dependence of the Gaussian 
component of the spin-spin relaxation rate T2G of a bi- 
layer system consisting of two different layers with doping 
xi = 0.16 and X2 = 0.24, respectively {5 = 0.4). For 
comparison, we show also the result for two bilayer systems 
with equivalent layers with doping xi = X2 = 0.16 and 
x\ = X2 = 0.24, respectively. The other parameters for all 
curves are: t' — —0.38t,t" = —0.06t. The inset shows the 
ratio Tia/TiG for the case of inequivalent layers. 

We start the discussion with the transverse relaxation 
time T2G- In Fig. ^ we show I/T^g as a function of 
temperature for both planes of the system with inequiva- 
lent layers in comparison with the corresponding data for 
two corresponding bilayer systems with equivalent layers, 
one with the same doping as the first layers and one with 
the same doping as the second layer of the first system. 
These results for the different layers are obtained with 
Eq. ^^by substituting with element (11) or (22) from 
Eq. ]10| . The filled (open) symbols in this figure are re- 
lated to the system with inequivalent (equivalent) layers. 
Furthermore, the squares represent the data for planes 
with a hole doping of a; = 0.16 and the bullets planes 
with doping x = 0.24. From Eq. |l3| it is seen that a sys- 
tem with rather strong magnetism and hence large values 
of x^^(ct) exhibits large relaxation rates 1/T2G- This ex- 
plains the differences between the layers with x = 0.16 
and X — 0.24, i.e. with different doping levels, whereby 
the heavily doped plane shows a smaller relaxation rate. 
The most striking result is that data for a plane with a 
given doping also depend on whether that plane is cou- 
pled with an equivalent one or with a more or less doped 
one. The heavily doped plane {x2 = 0.24) of the system 
with inequivalent layers shows stronger magnetic fiuctu- 
ations than the plane in the corresponding system with 



equivalent bilayers {xi = X2 = 0.24). Similarly, the mag- 
netism of the lower doped plane {xi = 0.16) is reduced 
with respect to the corresponding equivalent-layer sys- 
tem (xi — X2 — 0.16) due to the coupling to a stronger 
doped plane. Thus, the magnetic fiuctuations of the two 
inequivalent planes with different carrier concentration 
tend to be equalized by interplane coupling effects. A 
related effect has been detected also in Quantum-Monte- 
Carlo simulationSpOf coupled Hubbard planes carried out 
by Scalettar et al.Ej, in which the two planes have differ- 
ent chemical potentials and one plane is adjusted to half 
filling. In this case, the antifcrromagnetic susceptibility 
of the half-filled plane is reduced by the coupling with the 
doped plane. The authors explain this by the fact that 
processes in which holes hop from the doped layer into 
the half filled one are energetically more favorable than 
virtual hopping processes associated with the magnetic 
exchange J (x /U. 
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FIG. 12. Effective interaction V't,fr(Q,a; = 0) for a sys- 
tem with inequivalent layers {xi = 0.08 and X2 ~ 0.11) 
in comparison with the corresponding bilayer systems with 
xi = X2 ~ 0.08 and xi — X2 ~ 0.11. The upper panel shows 
the results for parameter set A and the lower for B. 

The effects discussed above are also visible in other 
quantities, like the effective interaction T4ff(q, w), which 
is proportional to the spin susceptibility in our approxi- 
mation. For this quantity, we can see that the stronger 
the antifcrromagnetic fluctuations are in the planes, the 
stronger is the tendency to equalization of the two planes. 
In Fig. ^(a) we show Vcsifi — {t^,t^),^ = 0) as a func- 
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tion of temperature for the parameter set A as measured 
in the two layers of the system with inequivalent layers. 
In the same figure, we also report for comparison the data 
for the two corresponding systems with equivalent layers 
with their doping adjusted to the one of each of the two 
layers of the first system. In Fig. |l^(b) we show the same 
comparison for the systems with parameter set B. The 
difference is striking. For parameter set A, VeS of each 
of the two layers in the system with inequivalent layers 
is considerably different from VeS in the corresponding 
system with equivalent layers and tend to be equalized 
for the two layers. On the other hand, for the less anti- 
ferromagnetic parameter set B, Veg calculated on a given 
layer of the system with inequivalent layers is essentially 
the same as the one calculated on the system with the 
same doping and equivalent layers. 
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FIG. 13. Intra- and inter-plane spin-spin relaxation rates 
I/T20 and I/T20, respectively. The inset shows the ratio be- 
tween the two relaxation rates {xi — 0.16, X2 = 0.24, S = 0.4). 

Qualitatively, the experimental situation for the sys- 
tem with inequivalent layers, Y2Ba4Cu70i5, appears to 
be described by pararaeter set A. Indeed, NMR exper- 
iments by Stern et al.H on I/T2G for the two layers of 
Y2Ba4Cu70i5 and for the two associated systems with 
equivalent layers, YBa2Cu307 and YBa2Cu408, show 
the same behavior as observed in Fig. |ll| and in Fig. 
Fig. |l^(a), if one identifies the behavior of I/T2G with 
the one of Vcs- Both experimental and theoretical re- 
sults show a strong increase of I/T2G as the temperature 
is lowered. However, we do not obtain the decrease of 
1/72G below Tsg ~ lOOK which is attributed to the open- 
ing of a gap in the spin excitation spectrum, since this 
region is probably unaccessible by our approximation. 
Another experimental observation is that the spin-lattice 
relaxation rate I/T2G has the same temperature depen- 
dence in the two planes of Y2Ba4Cu70i5. This has been 
deduced from the ratio R = {IIT:^'g)I{IIT^g) of 1/'^2G 
in the two planes, which turned out to be temperature 
independent and approximately R « 1.4-1.5. Since the 
124 plane in the coupled layer structure of Y2Ba4Cu70i5 
is the one with lower doping, the Cu02-layer from the 123 



block corresponds to the second plane in our theoretical 
study. The calculated values for R = T^'' /t!^'' are 
presented in the inset of Fig. |l^. It turns out that R is 
almost independent of the temperature T, in agreement 
with the experimental finding. 




A and various values of xi and 5. 



In fact, a systematic study shows that the ratio R is 
mainly controlled by the on-site energy 5 and thus by the 
difference of particle densities, 5n = rii — n2, between the 
two planes. On the other hand, the temperature depen- 
dence of R{T) is sensitively related to the doping of the 
first layer: for low doping xi, R{T) shows an upturn with 
decreasing T while it is essentially constant for large val- 
ues of Xi (see Fig. |lj). Thus, the experimental data for 
R{T) imply a rather large doping xi of the first layer 
in connection with a large on-site energy difference S or, 
equivalently, filling difference Sn. On the other hand, 
the doping cannot be too large because otherwise no ten- 
dency to equalization would be observable. An optimal 
choice for the doping levels in the coupled system turns 
out to be xi = 0.16 and X2 — 0.24, as shown in Fig. |ri|. 



For the same parameter choice, we compare the tem- 
perature dependence of the in-plane relaxation rate 
I/TIq with the inter-plane one I/T^q, which is calcu- 
lated using Eq. |^ The latter quantity has been mea- 
sured using NQR|-SEDOR experimentso, as suggested by 
Monien and RiceEl. The apparent feature in these exper- 
iments is that the inter-plane relaxation rate increases 
faster for decreasing temperature than the in-plane one, 
as seen from the temperature dependence of the ratio 
^SEdor(7') = T^ n/T^n- Our theoretical Calculations dis- 
played in Fig. |l^ clearly reproduce the qualitative behav- 
ior observed experimentally. However, we do not observe 
the saturation effect for very low temperatures as seen 
in experiments and which is most probably due to the 
opening of the spin gap. 



9 



0.25 

0.20 

0.15 

0.10 

0.05 

0.00 
0.00 



0.16/0.24 
0.16/0.24 
0.16/0.16 
0.24/0.24 



1.5 I- TjT^^/TjT^' 




-B 



0.05 



0.10 
T/t 



0.15 



0.20 





0.05 




0.04 


p 




H 


0.03 




0.02 




0.01 




0.00 




0.00 



0.16/0.24 
0.16/0.24 - 
0.16/0.16 ■ 
0.24/0.24 

0.05 



■-E3 

-0-- 



0.10 
T/t 



0.15 



0.20 



FIG. 15. Temperature dependence of the spin-lattice re- 
laxation rate l/{TiT) of a system with inequivalent layers 
compared to the corresponding systems with equivalent lay- 
ers and the same doping concentration, the inset shows the 
ratio of l/(TiT) for the two layers of the inequivalent-layer 
system (cf . Fig. |ll]) . 
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A similar study of the spin-lattice relaxation time Ti 
on Cu-sites shows less clear coupling effects between the 
two layers, both theoretically and experimentally. The 
experimental resultsH show almost no difference between 
the l/(TiT) data for the 123/124 layers in Y2Ba4Cu70i5 
and the corresponding layers in the YBa2Cu307 and 
YBa2Cu408 systems, respectively. Deviations can be 
seen only for rather low temperatures T < Tsg. Above 
Tsg, the ratio Ri{T) = {TiT^^^) / (TiT^^^) deduced from 
the experiments is again essentially temperature inde- 
pendent and is about 1.4. In Fig. [ij we show the cal- 
culated results for l/(TiT) obtained with Eq. ^ for the 
same parameters as in Fig. |ll|. Here, we observe that 
the coupling effects between the planes are not as strong 
as for the spin-spin relaxation rate 1/T2G- In fact, es- 
pecially for the plane with the larger doping {x — 0.24) 
\/{TiT) seems to deviate from the value of the system 
with equivalent layers only for rather small temperatures 
T/t K, 0.05. On the other hand, the plane with lower dop- 
ing {x = 0.16) shows deviations already for higher tem- 
peratures T/t « 0.08. Even more surprising is the result 
for the ratio i?i(T), shown in Fig. |ll|. Here, the plane 
with the higher doping shows a larger spin-lattice relax- 
ation rate \/{TiT), in contradiction with the simple ex- 
pectation that lower doping should result in stronger an- 
tiferromagnetic fluctuations. The expected value i?i > 1 
is restored only for low T % 0.04t in our calculation. The 
deviation from the experimental ratio i?i « 1.4 may be 
caused by the incorrect assumption that the hyperfine 
interaction constants are the same in both planes. 



FIG. 16. Temperature dependence of the spin-lattice re- 
laxation rates l/i^'^TiT) and of the spin-spin relaxation rates 
!/{}'' T2g) on oxygen sites. Like in Fig. ^ and Fig. ^ the 
results for a system with inequivalent layers are compared 
with the corresponding systems with equivalent layers and 
the same doping (the parameters are the same as in Fig. p^ . 

Finally, we have calculated the relaxation times ^^Ti 
and ^^T2G on oxygen sites. Like the corresponding re- 
laxation times on copper sites, these quantities are re- 
lated to the slope for w ^ of the imaginary part of 
the spin susceptibility, and by its real part, respectively. 
The appropriate form factor which enters Eqs.|ll|,|l3| and 
which is given by ^"^ F{q) — 2C^[l + 0.5[cos(q^)-t-cos(qj^)]] 
with two different constants for TiT and TIjg, suppresses 
scattering events transfering momentum q = (tt, tt) re- 
lated to magnetic fluctuations and instead favors those 
with q around (0,0). The spin-lattice and spin-spin re- 
laxation times on oxygen sites hence mainly probe the 
center of the Brillouin zone. We already know from the 
static spin susceptibility shown in Fig. ^ that the equal- 
ization effects are much weaker for q = (0, 0) as com- 
pared with (tt, tt). This weaker connection between the 
layers at q = (0, 0) is equally pronounced for the dynami- 
cal susceptibility and thus leads to the oxygen relaxation 
rates l/("TiT) and l/("r2G) shown in Fig. |l|. Here, 
we find that the layers of the system consisting of two 
inequivalent layers behave almost like the corresponding 
planes in the systems with equivalent layers and the two 
inequivalent layers are essentially disconnected. We thus 
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predict that, within a purely magnetic scenario, experi- 
mental measurements of the oxygen relaxation times in 
123 and 124-layers of Y2Ba4Cu70i5 should behave, as a 
function of temperature, like those of the corresponding 
layers in YBa2Cu307 and YBa2Cu408, respectively. In 
other words, if one would detect a different behavior of 
l/(^^TiT) in Y2Ba4Cu70i5 compared to the two parent 
compounds, it would be a strong indication for a non- 
magnetic coupling of the two layers. 



V. CONCLUSIONS 

In summary, we have studied a microscopic model con- 
sisting of two inequivalent Hubbard planes which are con- 
nected by an interlayer hopping t±^ . We have shown that 
magnetic and single-particle fluctuations of the two lay- 
ers are connected and tend to be equalized if the anti- 
ferromagnetic fluctuations within the layers are strong. 
If the antiferromagnetic correlation length is less than 
1-2 lattice spacings, which happens for high tempera- 
tures, large doping or bandstructure parameters with in- 
efficient magnetic coupling between the two Fermi sur- 
face sheeds, we find that the two inequivalent layers are 
disconnected and keep their individual properties. How- 
ever, once the antiferromagnetic correlations in the layer 
with smaller charge carrier concentration is sufficiently 
large, the single-particle excitations for momenta close 
to the Fermi surface and, in particular, around the hot- 
spots as well as the magnetic excitations of the two layers 
are strongly connected. The whole system reacts, de- 
spite its inhomogenious charge density, magnetically as 
a single system. In this case, the interlayer antiferro- 
magnetic susceptibility, measured by T2Q, increases for 
decreasing temperatures more strongly than the individ- 
ual inplane susceptibilities. These trends are in agree- 
rrient with the experimenmtal observation by Stern et 
an demonstrating that it is sufficient, for an understand- 
ing of the magnetic interlayer coupling, to use a single- 
particle interlayer hopping element. Furthermore, we ex- 
pect from our analysis of the oxygen NMR relaxation 
rate that, in distinction to the Cu-relaxation rates, the 
magnetic connection between the two layers will barely 
be visible, even for strongly underdoped systems. This 
phenomenon could be used to separate the small contri- 
butions due to antiferromagnetic correlations from the 
dominant, rather conventional, contribution of the q « 
dynamical spin susceptibility. Finally, in a futher step, 
based on the findings of this paper, it is of interest to 
investigate the behavior of two coupled, but inequivalent 
layers in the superconducting state (some work in this di- 
rection, although at much higher temperatures, has been 
carried out in Ref. [l9| ). This may reveal, why the mate- 
rial Y2Ba4Cu70i5 exhibits a higher Tc than both of the 
corresponding parent compounds it consists of. 
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